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ACCOMPLISHMENTS  UNDER  THE  CONTRACT 

A.  STATEMENT  OF  THE  PROBLEM  STUDIED 

A  three  year  program  of  research  was  proposed  to:  1)  Utilize  recently  developed  WSU  fast  line- 
scan  techniques  to  study  the  temperature  distribution  in  the  vicinity  of  the  crack  tips  of  propagating  brittle 
fracture  cracks  in  structural  polymer  composites;  2)  To  utilize  Wayne  State  University's  well  established 
box-car  video  thermal  wave  imaging  technique  for  the  purpose  of  studying  adhesive  bond  strength,  and  3) 
Exploit  the  concept  of  vector  lock-in  video  thermal  wave  imaging  for  mechanical  and  tlicrmal  response  of 
semiconductors. 

B .  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

1)  The  Thermal  conductivity  of  Isotopically  Modified  Single  Crystal  Diamond. 

A  tlicnnal  wave  materials  characteri/ation  project  in  diamond,  with  ARO  support,  has  led  to  a 
fundamental  understanding  of  that  material  as  a  heat  spreader.  Our  results  were  published  in  a  Physical 
Review  Letter;  "The  Thennal  conductivity  of  Isotopically  Modified  Single  Crystal  Diamond."  Lanhua 
Wei.  P.K.  Kuo.  R.L.  Thomas.  T.R.  Antliony  and  W.F.  Banhol/er.  Phys.  Rev.  Leders  70.  No.  24. 
pp.  3764-3767  (14June.  1993).  This  re.search  was  the  culmination  of  a  collalxiration  of  Wayne  State 
University  and  General  Electric  Corporate  R&D.  in  which  we  found  new  exi)crimental  results  on  the 
themial  conductivity  of  isotopically  enriched  '^C  diamond  crystals  at  low  temperatures.  To  our 
knowledge,  the  measured  value  for  a  99.9%  '^C  crystal  at  104K.  410W/cmK.  is  the  highest  measured 
thermal  conductivity  for  any  solid  above  liquid  nitrogen  temperature. 

Diamond  as  a  material  has  many  outstanding  physical  properties,  one  of  the  most  striking  of  which 
is  its  thermal  conductivity  (5  times  that  of  copper  at  room  temperature).  The  two  stable  isotopes  of  carbon, 
'^C  and  '^C,  have  natural  abundances  of  98.9%  and  1.1%,  respectively.  It  was  reported  by  the  authors  in 
1990  that,  by  enriching  the  isotopic  purity  of  '^C  from  the  natural  98.9%  to  99.9%,  the  thermal 
conductivity  of  diamond  crystals  at  room  temperature  is  enhanced  by  nearly  50%.  The  magnitude  of  the 
enhancement  is  surprising,  since  in  similar  systems  such  as  LiF  (*Li  and  ’Li)  or  He  (  He  and  He)  the 
isotopic  effect  is  only  1-2%  at  the  same  reduced  temperature.  The  unusual  enhancement  in  the  thermal 
conductivity  of  diamond  at  room  temperature  has  since  been  confirmed  by  researchers  from  other 
laboratories  around  the  country  using  a  variety  of  techniques.  The  Wayne  State  researchers  have  also 
reported  the  dependence  of  the  room  temperature  thermal  conductivity  on  content  over  a  range  from 
0.07%  to  99%  '^C. 

In  Fig.  1.  [from  the  Phys.  Rev.  Letter  cited  above],  we  summarize  the  results  of  fits  of  our  data  to 
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Fig  1  Thermal  conducUvity  of  natural  abundance  (1,1%  ’’C)  diamond  (lower  squares),  isotopicdly 
enriched  (0.1%  ’’C)  diamond  (upper  squares),  together  with  the  low  temperature  data  of  Slack  (circle) 
and  the  high  temperature  data  of  Olson  et  al.  (crosses).  The  solid  curves  are  the  results  of  fitting  the 
Callaway  theory  to  the  data,  using  the  same  set  of  fitting  parameters.  The  dashed  curves  arc  the  results  ot 
fitting  with  the  assumption  that  N-processes  dominate.  The  inset  shows  the  calculated  thermal 
conductivity  corresponding  to  1%,  0.1%,  and  0.001%  ‘’C-concentrations  according  to  the  Callaway 
theory. 


the  full  Callaway  model,  extended  to  take  into  account  the  differences  among  the  longitudinal  and 
transverse  phonon  modes.  In  the  paper,  and  as  illustrated  in  Fig.  1.  we  show  that  our  measured 
temperature  dependent  conductivities  for  the  isotopically  enriched  diamond  and  natural  abundance  diamond 
specimens  are  well  described  by  Callaway's  theoretical  model.  We  predict  that  the  thermal  conductivity  of 
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a  99.999%  C  diamond  crystal  should  exceed  2000  W/cmK  at  -80K.  It  is  of  practical  importance  that 
this  peak  is  predicted  to  occur  above  liquid  nitrogen  temperature. 

2)  Thermal  Wave  Imaging  of  the  Temperature  Distributions  around  Propagating 
Cracks  in  Polymers. 

The  research  project  was  carried  out  in  collaboration  with  Dr.  Henk  Van  Oenc  of  Ford  Scientific 
Research  Labs.  It  was  aimed  at  examining  the  spatial  and  temporal  distribution  of  energy  dissipation  in  the 
vicinity  of  brittle  fracture  cracks  in  polymer  materials.  The  technique  utilized  infrared  (IR)  cameras, 
operating  in  two  different  spectral  bands  (3pm  -  5pm  and  8pm  -  12pm)  of  the  IR,  and  witli  different 
image  frame  rates,  ranging  from  30  Hz  to  244  Hz.  A  table-top  tensile  tester  was  used  to  apply  tensile 
stress  to  notclied  polymer  test  coupons,  using  unifomi  drawing  speeds  which  arc  varied  up  to  100 
mm/minute  (see  Fig.  2).  During  the  tensile  test,  the  IR  emission  from  Uie  coupon  is  monitored  as  a 
function  of  time,  using  a  SantaBarbara  128  x  128  InSb  focal  plane  array  camera  (3pm  -  5pm  .spectral 
band/ up  to  244  Hz  frame  rate).  By  way  of  illustration,  in  Fig.  3  we  show  a  set  of  selected  frames  from 
such  a  tensile  test  on  a  coupon  of  rubber-modified  poly.styrene.  using  the  Santa  Barbara  camera  at  244  Hz. 
In  Fig.  4.  we  show  one  frame  Uiken  at  the  initiation  time  of  a  crack  in  a  rubber-modified  polystyrene  test 
sample,  with  a  spatial  resolution  of  about  40pm/pixel.  It  can  be  .seen  from  Fig.  4,  tliat  Ure  position  of 
greatest  IR  emission  in  the  vicinity  of  the  crack  tip  can  be  measured  quite  precisely  as  a  function  of 
position,  with  a  temporal  resolution  of  about  4  msec  (with  this  camera).  Fig.  5  shows  a  later  frame  from  a 
similar  sequence  to  that  shown  in  Fig.  2.  illustrating  a  "snapshot"  of  the  IR  emission  profile  in  Uie  vicinity 
of  the  propagating  crack.  Tlie  general  features  of  tire  experiment  confirm  tliat  a  material  with  a  positive 
thermal  expansion  coefficient  cools  when  it  is  deformed  elastically  in  tension.  In  some  of  our 
experiments,  we  utilize  the  detection  of  this  elastic  cooling  to  trigger  the  data  acquisition  which  follows.  In 
other  cases,  we  trigger  from  a  simultaneous  stress-strain  measurement,  beginning  to  record  images  at  the 
onset  of  plastic  deformation. 


-  3- 


ON  MECHANICAL  PROPERTIES  OF  MATERIALS  BY  THERMAL  WAVE  IMAGING 
R.L.  Thomas.  L.D.  Favro.  and  P.K.  Kuo,  Wayne  State  University 


Frame  #20  Frame  #30  Frame  #40 


Frame  #50  Frame  #60  Frame  #70 


Frame  #80  Frame  #90  Frame  #1 00 


Fig.  3  Sequence  of  frames  taken  during  the  initiation  and  propagation  of  a  crack  in  a  rubber-modified 
polystyrene  test  sample,  using  the  Santa  Barbara  camera  at  244. 
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Fig.  4  One  frame  taken  at  tlie  initiation  time  of  a  crack  in  a  rubber-modified  polystyrene  test  sample, 
using  the  Santa  Barbara  camera  at  244  Hz,  and  with  a  spatial  resolution  of  about  40pm/pixel. 
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Fig.  5  One  frame  taken  during  the  propagation  time  of  a  crack  in  a  rubber-modified  polystyrene  test 
sample,  using  the  Santa  Barbara  camera  at  244  Hz.  with  a  spatial  resolution  of  about  120pm/pixel. 

We  assume  that  plastic  yielding  occurs  mostly  near  the  crack  tip,  so  that  the  location  of  the  crack  tip 
is  revealed  by  the  largest  observed  temperature  in  the  image  field.  A  typical  experiment  yields  about  a 
hundred  consecutive  images.  Hence,  a  rather  dense  set  of  data  can  be  analyzed,  both  for  the  initial 
yielding  and  crack  intiation,  and  for  the  subsequent  "catastrophic"  propagation.  An  example  of  the  raw 
data  for  three  samples  of  high-impact  polystyrene,  taken  at  three  different  pulling  speeds,  is  shown  in 
Fig.  6. 
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Crosshead  speed:  a  -  40  mm/min 

o  -  60  mm/min 

Fig.  6  Plots  of  the  crack  length  as  a  function  of  time,  for  three  different  pulling  speeds,  using  the 
apparatus  of  Fig.  2.  The  crack  length  measurements  were  taken  from  the  positions  of  the  biggest 
temperature  pixels  in  the  sequence  of  stored  images. 

Together  with  Dr.  Henk  van  Oene  of  the  Polymer  Science  Department  of  the  Ford  Scientific 
Research  Laboratory,  we  have  recently  analyzed  these  crack  initiation  and  propagation  results  in  terms  of 
the  theory  of  Levy  and  Herrmann'.  For  an  arbitrary  load,  P ,  one  can  write 


'  C.  Levy  and  G.  Herrmann,  Eng.  Fract.  Mech.  17,  125  (1983). 
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Gda:^ 


P^dC 


(1) 


where  G  is  the  energy  release  rate  per  unit  crack  length;  a  =  1  +  b.  where  1  is  the  crack  length,  and  b  is  the 
notch  depth.  If  we  let  C  be  the  compliance  of  the  beam,  it  can  be  written  as 


P  2a 

E  .  (2) 


where  E  is  the  tensile  modulus  and  a  is  a  function  of  4  =  a/h.  and  h  is  the  width  of  the  beam.  According  to 
this  theory,  for  jt/li<0.6,  the  function,  a.  is  given  by  the  tenth-order  polynomial 


a  =  1.98  4-  -  0.54  4’  +  18.65  4^  -  33.7  4'  +  99.26  4‘  -  21 1.9  4’  -H136.84  4*  -460.77  4''  -289.98  4'°  .  (3) 

Taking  into  account  that  the  energy  release  rate  can  al.so  be  expressed  in  tenns  of  the  fracture  toughne.ss.  K, 


(4) 


one  .sees  from  Eq.  ( 1 )  that 


2  da 
"  h  d4 


(5) 


The  testing  is  carried  out  at  a  constant  rate  of  extension.  For  a  given  extension,  5,  we  have 

5  =  CP. 


(6) 


Hence, 


dt 


pdP  p^l 

da  da  dt 


(7) 


Since  the  rate  of  extension  is  held  constant  during  the  experiment,  on  may  derive  an  expression  for  E; 
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a  dP  P  da  da/dt 
Of  dioJOf)  ^  h  d5/dt 


(8) 


A  "dynamic"  energy  release  rate  can  now  be  defined  as 


G 


dynamic 


(9) 


where  is  calculated  from  Eq.  (5)  and  E  is  calculated  from  Eq.  (8).  once  Uic  crack  velocity,  da/dt.  is 
known.  Even  though  the  crack  velocity  could  be  calculated  from  the  data,  as  incremental  growth  divided 
by  the  frame  time,  this  procedure  will  yield  unsatisfactory  results,  since,  intitially,  the  the  crack  velocity  is 
very  slow  and  the  data  set  is  sufficiently  den.se  that  the  incremental  crack  advance  contitutes  a  distance  of 
O.l,  or  2  pixels.  Instead,  tlic  crack  tip  location,  as  revealed  by  the  pixel  with  the  largest  incrca.se  in 
temperature,  is  fitted  to  an  expression  for  the  crack  velocity  whicli  lias  licen  proposed  by  M.L.  Williams’ 

^  =  cxiXAt) 


from  which  ihc  crack  velocity  can  be  calculated  as: 

'  \/b  .  (11) 

From  a  plot  of  the  left  side  of  Eq.  (10)  as  a  function  of  time,  one  can  determine  the  value  of  the  contant  A 
at  different  pulling  speeds.  An  example  of  plots  such  as  this  is  shown  in  Fig.  7. 


^  M.L.  Williams.  "The  Fracture  of  a  Viscoelastic  Material"  in  Fracture  of  Metal  Solids,  Metallurgical  Society 
Conferences,  Vol.  20.  D.C.  Druckcr  and  J.J.  Gilman  (Eds  ).  Interscience  Publishers.  New  York  (1963). 
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Fig.  7  Plots  of  tlic  log  of  the  left  side  of  Eq.  (10)  as  a  function  of  time.  These  plots  can  be  used  to 
determine  the  constant  A  in  Eqs.  (10)  and  (11). 


The  values  of  the  contant  A  determined  from  Fig.  7  can  tlicn  be  inserted  in  Eq.  ( 1 1 )  to  obtain  a  scries  of 
plots  of  the  crack  velocity  as  a  function  of  crack  length.  Tlic  universality  of  Eq.  (10)  can  be  seen  in  Fig.  8, 
which  plots  1/b  versus  reduced  time.  At.  It  seen  from  Fig.  8  tliat  all  of  tlie  data  for  1/b  >  1  fall  on  such  a 
universal  curve.  Plots  of  the  crack  velocity  versus  crack  length  from  Eq.  (11),  using  these  values  of  the 
constants  A,  are  shown  in  Fig.  9. 
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Fig.  8  Plot.s  of  1/b  a.s  a  function  of  rcducce/ time.  At.  The  fact  that  all  of  the  plots  fit  on  a  single  curve 
demonstrates  the  universality  of  lEq.  (10). 
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Fig.  9  Plots  of  the  crack  velocity  as  a  function  of  crack  length,  using  the  values  of  tlie  constant  A 
determined  from  Fig.  7  and  Eq.  (10). 

In  Fig.  10,  we  present  plots  of  the  maximum  temperature  at  the  crack  tip,  as  a  function  of  time,  for 
five  different  pulling  speeds,  using  the  apparatus  of  Fig.  2.  The  temperature  measurements  were  taken 
from  the  values  of  the  highest  temperature  pixels  in  the  sequence  of  stored  images.  It  is  interesting  to  note 
that  the  temperature  rise  tends  to  saturate  after  an  initial  rapid  rise. 
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Fig.  10  Plot.s  of  (lie  maximum  temperature  at  the  crack  tip.  as  a  function  of  time,  for  five  different 
pulling  speeds,  using  tlic  apparatus  of  Fig.  2.  The  temperature  measurements  were  taken  from  the  values 
of  the  highest  temperature  pixels  in  the  .sequence  of  stored  images. 


In  addition  to  the  experiments  on  single-notched  coupons,  described  above,  we  have  also  imaged 
some  double-notched  coupons  to  study  the  interaction  between  two  propagating  cracks,  and  single- 
notched  coupons  with  defects  to  study  tlic  interaction  of  propagating  cracks  witli  defects.  In  Fig.  1 1 ,  we 
show  a  sequence  of  images  for  a  double-notched  coupon  of  high-impact  polystyrene,  which  show 
adiabatic  cooling  during  elastic  deformation,  followed  by  crack  tip  initiation  and  plastic  deformation, 
propagation  of  the  two  cracks,  and  the  interaction  of  the  two  crack  tips  prior  to  fracture.  The  fact  that  the 
two  crack  tips  are  seen  to  "repel"  one  another  is  quite  characteristic,  and  relates  to  the  stress  relief  in  the 
plastic  deformation  zone  surrounding  each  crack. 
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Fig.  1 1  Sequence  of  images  for  a  double-notched  coupon  of  higlyimpact  polystyrene,  which  show 
adiabatic  cooling  during  elastic  deformation,  followed  by  crack  tip  initiation  and  plastic  deformation, 
propagation  of  the  two  cracks,  and  the  interaction  of  the  two  crack  lips  prior  to  fracture.  Tire  fact  tliat  tlie 
two  crack  tips  are  seen  to  "repel"  one  anotlier  is  quite  characteristic,  and  relates  to  tlie  stress  relief  in  lire 
plastic  deformation  zone  surrounding  each  crack. 


In  Fig.  12,  we  show  a  sequence  of  images  for  a  single-notched  coupon  of  polycarbonate,  containing  a 
defect  in  the  form  of  a  bubble.  The  high  elastic  stress  region,  both  at  the  tip  of  the  notch,  and  around  the 
bubble  ate  seen  ftrst  to  cool,  and  then  to  heat  as  the  plastic  deformation  sets  in.  While  the  plastic 
deformation  region  passes  through  the  bubble,  the  crack  docs  not.  This  behavior  is  consistent  with  that 
seen  for  the  double-notched  example  of  Fig.  1 1,  showing  that  cracks  avoid  regions  in  which  the  stress 
has  been  relieved. 


-  15- 


ON  MECHANICAL  PROPERTIES  OF  MATERIALS  BY  THERMAL  WAVE  IMAGING 
R.L  Thomas.  L.D.  Favro.  and  P.K.  Kuo.  Wayne  Stole  University 


Locnl  Hcniing 


ts40l  ms 


t 

k 

Crack<Flaw 

Inferi^ctton 


Fig.  12  Sequence  of  images  for  a  single-notclicd  coupon  of  polycarbonate,  containing  a  defect  in  die 
fonn  of  a  bubble.  Ilie  high  elastic  stress  region,  both  at  the  tip  of  the  notch,  and  around  die  bubble  arc 
.seen  first  to  cool,  and  then  to  heat  as  the  plastic  deformation  sets  in.  While  the  pUislic  defomiation  region 
passes  through  the  bubble,  the  crack  docs  not.  This  behavior  is  consistent  with  diat  seen  for  die  double- 
notched  example  of  Fig.  1 1.  showing  that  cracks  "avoid"  regions  in  which  the  stress  has  been  relieved. 

3)  Thermal  Wave  Imaging  of  the  Fracture  of  adhesive  joints. 

Tlie  third  topic  for  which  we  present  important  results  is  that  of  fracture  in  adhesive  bonds.  In 
Fig.  1 3,  we  show  four  images  of  a  lap  joint  between  two  metal  plates,  viewed  from  the  side  of  the  joint, 
as  the  joint  undergoes  tensile  stress  to  failure.  It  can  be  seen  that  cracks  intiate  at  both  ends  of  the  joint 
and  propagate  towards  the  center,  with  a  large  energy  release  as  the  final  fracture  of  the  joint.  This  work 
was  carried  out  in  collaboration  with  Mr.  A.C.  Ramamurthy  of  Ford.  A  somewhat  different  adhesive 
fracture  experiment  was  carried  out  in  collaboration  with  Dr.  Ray  Dickie,  also  of  Ford.  Some  of  this 
more  rccent  work  has  been  carried  out  with  a  faster  frame  rate  (1  kHz)  IR  Camera.  .  An  expenmental 
arrangement  for  imaging  such  fracture  events  is  shown  in  Fig.  14. 
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Fig.  13  Sequence  of  four  iinagc.s  of  a  lap  joint  iKlwccn  two  metal  plates,  viewed  from  the  side  ol  the 
Joint,  as  the  Joint  undcrgoc.s  tensile  stress  to  failure.  It  can  be  seen  that  cracks  intiatc  at  both  ends  ol  Uic 
joint’and  propagate  towards  the  center,  with  a  large  energy  relea.se  as  the  final  fracture  of  the  Joint. 
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Experimental  Setup  for  Studying  Propagating 
Cracks  in  Adhesively  Bonded  Metal  Plates 


Fig.  14  lExpcrimcnlal  .setup  for  studying  projxigating  cracks  in  adhesively  bonded  metal  plates.  A  fast 

IR  camera  images  the  joint  from  the  edge,  focussed  just  above  the  tip  of  tlie  wedge  where  the  crack  in  tltc 
adhesive  initiates  and  propagates. 
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Fig  15  Two  images,  taken  10  msec  apart,  of  the  propagation  of  a  crack  in  a  metal/metal  adliesivc 
bond,  using  the  experimental  setup  shown  in  Fig.  14.  The  total  time  from  initiation  of  the  crack  to  fracture 

was  18  msec. 


We  have  recently  acquired,  through  an  NSF  Instrumentation  Grant,  combined  with  University  and 
Industry  cost-sharing,  a  512  x  512  InSb  focal  plane  array  camera,  with  unique  windowing  and  high 
frame  rate  capability.  Tlie  frame  rate  for  the  full  aaay  is  up  to  72  Hz.  In  addition,  tlierc  are  eight  mouse- 
selecuible  windows,  with  higher  fnune  rates,  up  to  34  kHz.  Furthermore,  the  array  operates  in  a 
"snapshot"  mode  -  a  feature  of  importance  for  simplifying  the  analysis  of  the  time-dependence  for  various 
features  in  the  image  field.  We  have  also  recently  acquired  a  new  tabletop  tensile  tester,  with  a 
temperature-controlled  sample  cell  and  IR-transparent  window,  which  will  permit  us  to  extend  the 
measurements  described  to  temperatures  otlier  than  room  temperature.  We  are  also  using  this  camera  to 
continue  our  lock-in  imaging  experiments  on  semiconductors. 

Wayne  State  University  patents  and  other  intellectual  property  in  the  area  of  Thermal  Wave  Imaging 
(supported  in  part  by  ARC),  are  currently  undergoing  technology  transfer  by  a  start-up  company; 
Thermal  Wave  Imaging.  Inc.  (TWI),  18899  West  Twelve  Mile  Rd..  Lathrup  Village.  MI  48076  [Phone: 
(810)  569-4960;  FAX:  (810)  569-4252].  TWI  was  formed  to  license  from  WSU  and  commercialize 
patents  held  by  WSU  in  the  area  of  thermal  wave  imaging.  This  technology  transfer  activity  has  led  to  the 
EchoTherm'image  acquisition  hardware  and  software  (for  IBM-PC  compatible  computers),  which  it  is 
currently  offering  for  sale.  TWI  is  akso  engaged  in  consulting  activities  for  companies  interested  in 
applying  this  existing  thermal  wave  imaging  technology  to  their  proprietary  NDE  problems. 
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